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Abstract. We present a constructive approach for approximating the conformal map (uniformization) of a 
polyhedral surface to a canonical domain in the plane. The main tool is a characterization of convex spaces 
of quasiconformal simplicial maps and their approximation properties. As far as we are aware, this is the 
first algorithm proved to approximate the uniformization of general polyhedral surfaces. 



1. Introduction 

A polyhedral surface S is defined by stitching planar polygons along congruent edges. A polyhedral surface 
can be endowed with a conformal structure making it a Riemann surface |23l [5] . The celebrated uniformiza- 
tion theory [31 [TU] then implies the existence of a conformal map : S — > T between S C M. d (d is typically 
3 in applications) and a topologically equivalent domain in the plane TcC. The main focus of this paper 
is building constructive approximations of this map for polyhedral surfaces with generally shaped polygonal 
faces. Topologically, we restrict our attention to disk-type surfaces. To date, we are not aware of any other 
existing algorithm that is proved to converge in the limit to the uniformization map for general polyhedral 
surfaces (i.e., with polygonal faces of arbitrary shape). 

The problem of approximating planar conformal mappings is considered well-understood and there is a 
wealth of methods that produce approximations to conformal mappings between planar domains |19[ I18[ [9] . 
Circle-packing [23 , imitates conformal mappings by replacing infinitesimal circles with finite one, and was 
shown to converge to conformal mappings as the circles are refined [201 [Ml 12] • However, it seems no full 
generalization to polyhedral surfaces with generally shaped faces exists [33J . Other constructions of "discrete 
uniformization" exist in the field of Discrete Differential Geometry (DDG) [5] where the focus is building a 
consistent and rich discrete theory. An example for such construction that can be used to compute discrete 
uniformization is by Springborn and coauthors 22, 4J. Discrete Ricci flow is another example |13j . Other 
popular constructions can be found in |16Ulll[2"l"lll4l [T] however no proof of convergence to the uniformization 
map is provided for any of these methods so-far. 

In this paper we construct algorithms for approximating the uniformization map with guarantees (i.e. with 
convergence proof). Without loosing any generality we can subdivide each polygon in S into triangles and 
henceforth assume we have a triangulation S — (V, E, F), where V = {vi} C M. d , is the set of vertices; 
E = {e^}, the set of edges; and F = {fj}, the set of oriented planar faces (triangles). 

We will construct successive approximations to \& over a series of regular subdivisions 
of the surface S = S Q -< S 1 -< ... -< S q . By regular subdivision we mean that at each 
level we cut every face into four similar faces by connecting the mid-edges points, see 
inset figure, and Figure[l](top row). Our approximations will be simplicial mappings 
$ 9 £ J- s<1 of S q , namely piecewise-afnne (over faces) and continuous mappings into 
the complex plane C, the collection of such mappings over S q will be denoted by J 75 . <I> 9 will converge 
locally uniformly to ^. We will restrict our attention here to topological disks, and take as a canonical 
domain T the equilateral A(ii, £2, £3) defined by its three corners tt = e i (*-i)2t/3^ = 12) 3. To set a unique 
target uniformization map we will mark three distinct (positively oriented) boundary vertices v±, V3 G V 
that will be mapped to the corners t\, t%, £3 (respectively). This fixes all the degrees of freedoms of the map. 
Figure [T] shows examples of simplicial approximations $ 9 to ^ : S — > T for a series of three refinements of 
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Figure 1. Approximation of the uniformization map with simplicial maps on subdivided 
versions of the origianl polyhedral surface S. 

a particular polyhedral surface: the middle row shows planar checkerboard texture mapped by the inverse 
of the simplicial maps $ 9 to visualize the "conformality" of the approximations. The bottom row shows the 
homeomorphic image of S q under <I> 9 onto the equilateral domain. Bright-red color indicates high conformal 
distortion, while grey indicates low conformal distortion. Note that the approximations are improving as the 
mesh is rchncd. The main result on which we build upon when developing the algorithms in this paper is: 

Theorem 1.1. Let ^ : S — > T be the uniformization map of a disk-type polyhedral surface to the equilateral 
T taking three prescribed boundary vertices of S to the corners of T '■ Let S q be the q th -level subdivided 
version of S. 

If, for an arbitrary but fixed e > 0, the argument of "J' is known up to an error o/±(^ — e), then one can 
construct a series of non-empty convex spaces U q = {<& q } C J-" 5 ' of simplicial maps of S q such that: 

(1) Every map $ 9 € U q is K-quasiconformal (QC) homeomorphism that maps S onto T, with some 
constant K independent of q. 

(2) Every series {& q } q >Q> where $ 9 £ U q , converges locally uniformly to the uniformization map ^ . 
That is, $ 9 o ty^ 1 ; 7~ — > T converges uniformly in any compact subset of Interior{T) to the identity 
map Id -T — > T '. 

Let us clarify the assumption "the argument of 'i'' is known up-to an error of ±(? — e)". What we mean by 
that is that at every point p G S, we can choose (arbitrary) chart (ifi, f2j), ifi : fli C S — > C, where p € $1^, 
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z — <Pi(p), and that we can point an angle r in the range 

(arg [(* o ^)'(z)] - (| - e) , arg [(* o ^)'{z)} + (| - e)) , 

for e > arbitrary small but fixed. Intuitively, if we know in which "half of S 1 (S 1 denotes the unit circle) 
the argument of the derivative of the map we are looking for resides in, then we can approximate ^/ via a 
convex program, namely looking for an element in a known convex subset of the simplicial maps U q C J- sq . 

This theorem will be proved in several parts: first, in Section [3] we will show that J- s<1 contains at-least one 
simplicial map that is: 1) quasiconformal, namely a homeomorphism with bounded conformal distortion, 2) 
its conformal distortion is converging to 1 with almost linear rate. Second, in Section |4j we will show that 
any series of quasiconformal simplcial maps $ 9 that their conformal distortion converge to 1 converges to the 
uniformization map VP as described in Theorem |1.1| In Scction[5]we will study the space of AT-quasiconformal 
simplicial maps J^f? c J- sq and characterize a collection of convex subsets { Ffr T } , r = (ti , .., rips 1 ) , Tj € S 1 . 
We will further show that given an ^ — e approximation of the argument of 'J', one can single out one of 
these convex subsets that will be non-empty. This will finish the proof of Theorem 1.1 

Building upon Theorem |1.1| we will suggest two algorithms. The first algorithm is exhaustive but theoretically 
fully justified: it will test many candidates r, one of which is guaranteed to lead to a non-empty space J~^ T - 
The number of candidates will be rather large but shown to be independent of q (i.e., exponential in the 
number of original faces |F| — |F°|). Due to the large number of candidates r this algorithm will mainly 
have theoretical importance but only limited practical applicability. Nevertheless, as far as we are aware this 
is the first algorithm to approximate the uniformization map for general faced polyhedral surfaces. 

The second algorithm will be greedy in nature: it will start with some arbitrary convex space J-f? T0 , and 
will iteratively move to other convex spaces J-f c q Tl ,J-f, q T2 ,... When terminating successfully it guarantees 
an approximation to the uniformization map, enjoying all the properties of the first algorithm (and Theorem 
[Li] ) . However, the drawback here is that we do not have a proof that it will always ends up with a non-empty 
space (i.e., terminate successfully). Nevertheless, it works well in practice and is much more computationally 
efficient than the (first) exhaustive algorithm. 

We start by defining the smooth and discrete conformal structures of S, set notations and a few preliminary 
lemmas. 

2. Smooth and discrete conformal structures 

As noted above, a polyhedral surface S admits a smooth (classical) conformal structure. We start by defining 
it. 

A conformal structure is defined by providing a conformal 
atlas, that is, a set of analytic coordinate charts [10]. We 
will define such an atlas already customized to our later 
constructions. We fix a constant < C < 1/4, and for 
each vertex Vi € V, denote by Mi the set of indices of 
1-ring neighbors to vertex i. 

We distinguish three types of vertices in our surface S: 
interior, boundary, and corner. Interior vertices are ver- 
tices in the interior of S, boundary are on the boundary 
but not one of the three corners i>i, ^2, V3. 

Interior charts. Let Vi S V be an interior vertex. Set 
Qi = 2n. Set the domain fij C S to be the interior 
of the convex hull (over the surface) of the points C,Vi + 
(1 — ()vi>, where i' 6 A/i (see inset, where £li for one 
particular vertex is colored in green) . We define the chart 
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(fi : Qi —¥ C by first rigidly unfolding each triangle in 
Qi onto the plane, taking Vi to the origin, and second, 

composing each (now planar) triangle with the map z i— > z li , 7, = ^ where 0j = y\ 0j j is the sum of angles 
at vertex Vi, where 6{j is the angle of face j (adjacent to vertex i) at vertex i. This composition is made 
(possibly by incorporating rigid transformations in the plane) such that ipi is globally continuous, see [25115] 
for more details and the inset for an illustration. Let us denote Aj = ^j(f2j) C C. 

Boundary charts. Let Vi £ V be a boundary vertex (i ^ 1, 2, 3). Set Qi = tt. The domain fl; is defined as 
the interior of the convex hull of Uj and C v i + (1 — C) w i': G and <^>j : fl^ — >■ C is defined via mapping the 
neighborhood congruently to the plane as before, composing it with z v- > z 7i , where 7» = §f . 

Corner charts. For the three corner vertices Ui,U2> w 3 we se ^ ®i = f an d define fij, i — 1,2,3, similar to 
the boundary vertices' charts with the exception of using the mapping z 1— > z 7i , where 7i = §N 

All the transition maps o i/?" 1 are conformal and therefore define a conformal structure over S : 
Lemma 2.1. The transition maps (p iA i — ipi o ipT; 1 are holomorphic. 

Proof. The holomorphy of the transition maps ipij/ — fi o ip^, 1 can be understood (see e.g. ,[5]) from the 
fact that these maps are compositions of similarities z 1— > az + S and the analytic maps (note that we avoid 
the origin) z 1— > z f>7r / 8i , ?9 = 1/3, 1, 2. The holomorphy across edges can be verified using standard extension 
theorems of conformal maps. □ 

Now that we have a (smooth, classical) conformal structure, the notion of conformal mappings from S to the 
plane C is well-defined; a map ^> : S — > C is conformal if for every chart (f2,-,y>j), the map ipi = 'J/ o cp7 
defined over Aj = </?,-(fij) is conformal in the classical sense. 

The main object of this paper is to approximate the conformal mapping \& that maps the polyhedral surface 
bijectively to the equilateral triangle T. The existence and uniqueness of such a mapping is set by the 
uniformization theorem [31 [TO] : 

Theorem 2.2. There exists a homeomorphism ^> : S — > T that is conformal in the interior of S , ^> : 
Interior{S) — > Interior(T^\ The map ^ is uniquely set once required to take Vi,V2, W3 to the corners h,t2, £3 
ofT. 

We will approximate \& by constructing quasiconformal simplicial mappings $ 9 : S q — > T from subdivided 
versions of S to T. 

The subdivided triangulations S q = (V q , E 9 , F q ) , q = 0, 1, 2, ... of S are 
constructed by the standard 1 — 4 subdivision rule in each triangle, where 
(S° = S. For example, in the inset the red mesh shows S , and we show 
in purple all the faces fy € F 2 n /j inside one face fj € F of S. 

We now turn to define the discrete conformal structure over S q . A discrete 
conformal structure is basically assigning angles to the corners of each 
triangle such that the angle sum of each triangle is 7r. Equivalently, we 
can embed each triangle in the Euclidean plane and think of it up-to a 
similarity transformation. A consistent discrete conformal structure will approximate the smooth one as 
the surface is refined. A simple way to do it is by mapping each triangle to the plane with the charts ifi 
and taking its image Euclidean triangle to define its conformal structure. Obviously, using different atlas, 
or assignments of triangles to charts will lead to different discrete conformal structure, but as we will show, 
at the limit it won't matter. So we will use the Atlas defined above, and describe an arbitrary assignment 
of faces to charts in F 9 . 

We start by associating, for q = 2, each triangle fj € F 2 to some chart {Qi, tpi}, i = ij : by making sure that 
Closure (fj) C Hi- 




Interior(Q) denotes the interior of the set Q. 
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For example, we can use the assignment rule shown in the inset figure (each color 
indicates association of that colored face to a different vertex of the triangle and 
hence to a difference chart). 

For any fj £ F q 7 q > 2 we associate the chart based on the triangle's (unique) 
ancestor in level q = 2. We define the discrete conformal structure for each face 

fj £ F q ,q > 2 by mapping the face's vertices Vj 1 ,Vj 2 ,Vj 3 to the plane, that is zg — (pi(vj e ),£ = 1,2,3, 
{fli,<fi} the associated chart i = ij, and using the Euclidean triangle A(zi, z 2 , z 3 ) to define the discrete 
conformal structure of fj. The inset below shows in blue (bottom-left) the triangle A(z\, z 2 , 23). 

Once we defined the discrete conformal structure we 
can define the discrete conformal distortion of a sim- 
plicial map $ 9 : S q -> C. We look per face fj £ F q , 
q > 2, and denote by Aj the affine map mapping the 
triangle A(zi, z 2 , z^) (remember that zg = <Pi{Vj l )) 
to the triangle A{<S> q (v h ), <f> q (v 32 ), $ q (v j3 )). We will 
also refer to Aj as $ 9 |/. expressed "in the local co- 
ordinate chart". Then, the discrete conformal dis- 
tortion D d ($« |/.) of <S> q |/. is defined by 

D d (<F \ fi )=B(A 3 ), 
where T)(Aj) denotes the standard conformal distor- 
tion of the planar affine map Aj, namely, T)(Aj) = 
— , the ratio of the larger to smaller singular values 
> Oj > 1 of the linear part of the planar affine 
map Aj. See the inset figure for an illustration. 

The discrete conformal distortion of the full simpli- 
cial map $ 9 is accordingly defined by 

(2.1) 




D d (<F) = max D d ($« \ fj ). 

For later use, we denote by T) d n(§ q ), D, C 5, the discrete conformal distortion of the map $ 9 restricted only 
to triangles contained in the set f2, that is 



(2.2) 



D d a ($ 9 ) 



max D d L4 7 -) 



Let us denote by J 75 ' the space of simplicial maps (i.e. , continuous piecewise affine) mapping the triangulation 
S q to the plane, satisfying the boundary constraints of mapping dS q (dS q denotes the boundary of the 
polyhedral surface S q ) bijectively to dT and taking the corner vertices wi,W2,«3 to the triangle's corners 
ti,t 2 ,t 3 . Further define the subset C J- s " , K = (Ki,K 2 , Ki F qi) to include only orientation-preserving 
homeomorphisms $ 9 £ J-^ such that D d ($ 9 \f.) < Kj. Denote |K| = maxifj, and by J-f? we will mean 
^\k.\<kF^ , that is, simplicial maps with maximum bound K on its conformal distortion per face. 

There are two technical issues to be taken care-of for later constructions. First, let us highlight a small 
technical property of the charts and the association rule of faces we have defined that will be used in 
Section [3] when we approximate the uniformization map with simplicial map: we show that there exists some 
positive gap between the triangles and the boundary of the charts' they are associated with. Denote the 
disk V(z, r) — {w \ \w — z\ < ?'}. 

Lemma 2.3. There exists some constant Eq > 0, such that 
for all fj £ .F 9 , q > 2, associated with chart {17^, ifi}, i = ij 



Proof. It is enough to prove the inclusion for q — 2. Take arbitrary face fj £ F = F , and one of its vertices 
Vi £ V = V°. Since when defining fii we took < ( < 1/4 we see by the association procedure that every 
triangle /,/ £ F q is far from Dili at least by some constant Ej > 0. When applying tp i: since it has bounded 
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derivatives (when considering it as map from the Euclidean faces to the complex plane) away from the vertex 
Vi, we have some other e'a > bounding this distance from below. Taking the minimum s'j over the (finite) 
set of faces in F and their vertices we finish the argument. □ 

A second issue is that we need to make sure that approximating the power maps h : z i— > z 7 , 7 > 
with simplicial maps (via sampling at the vertices and extending linearly) are if-quasiconformal with some 
universal bound on their distortion K > 1, independent of subdivision level q. For an example of a simplicial 
power map see Figure [2] 




Figure 2. The simplicial map h q ,q = 2 sampled from h(z) 



We conjecture the following: 

Conjecture 2.4. Let T = A(^,0,ry) be a triangle and denote the angle 8 — /(£, 0, 77). Further let T q be 
the q th level of regular 1-4 subdivision of T . Denote by h(z) = z 1 the power map, and assume that 7$ < ir. 
Then, the simplicial maps h q defined by sampling h{z) over the vertices of T q and extending by linearity are 
homeomorphisms that satisfy D(h q ) < K for some K > 1 independent of q. 

For our needs it is enough to prove the following, slightly weaker, result: 

Lemma 2.5. Let T = A(£,0, ry) be an isosceles triangle = \rj\) and denote the angle 6 = Z(£,0,7y). 
Further let T q be the q th level of regular 1-4 subdivision of T . Denote by h(z) = z 1 the power map, and 
assume that \"f\B < ^, and that 6 < 60.4°. Then, the simplicial maps h q defined by sampling h(z) over 
the vertices of T q and extending by linearity are homeomorphisms that satisfy D(h q ) < K for some K > 1 
independent of q. 



The proof for this lemma is rather technical and therefore deferred to Appendix \K\ 

Building upon this lemma, we can (without loosing generality) 
subdivide each of the original triangles fj e F of the polyhe- 
dral surface S, so that the two conditions \j~\9ij < £, and 
9ij < 60.4°, where 7 = are satisfied. We can also guaran- 
tee that every triangle that touches one of the original vertices 
V is isosceles (all the other triangles have all flat vertex angle^] 
and so their charts are rigid congruencies and hence quasicon- 
formal). Such subdivision is shown in the inset. Note that this 
subdivision produces a conformally equivalent polyhedral sur- 
face to S with identity as the conformal equivalence, and there- 
fore approximating the uniformization for this new polyhedral 
surface is equivalent to approximating the uniformization of 
the original polyhedral surface S. 




By flat vertex angles we mean the angle sum around a vertex is 2ir. 
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3. Feasibility 

We are now ready to prove our first result. Namely, that a subdivided triangulation S q , q > 2, can be 
mapped to T with a simplicial homeomorphism $ 9 £ where the discrete conformal distortion of the 
triangles is controlled. 

Theorem 3.1. There exist simplicial maps $ 9 £ Ti 9 , for K= (K u ...,K\g,\), Kj = 1 + 0(2-"*«), where 
Ki — min{7i,l} and face fj £ F 1 is associated with chart i = ij. That is, <I> 9 : S — > T are orientation 
preserving homeomorphisms that satisfy for every face fj £ F q 

D d (® q I/,) = 1 + C(2" K ' 9 ). 
In particular, $ 9 £ Ffc > where K q = 1 + (D(2~ Kq ), and k = min^ {ni}. 

Furthermore, for every face fj £ F 1 , min^gz jargc^A, — arg o tp^ ) r (z)~\ + k2jr\ = O (2~ qK ), where Aj 
is $ 9 |fj in local coordinates, z is the centroid of the triangle A(zi, zi, z^), zg — fi{vj t ), £ = 1,2,3, and Vj e 
are the vertices of the face fj . 

Lastly, for any fixed domain l!c5 satisfying Closure(nf\n V=$, D d n ($ q ) = 1 + 0(2~^). 

We will need an auxiliary Lemma regarding approximation of conformal mappings with simplicial maps: 

Lemma 3.2. Let f : 0(0, r) — > 0(0, 1) be a conformal map such that |/'(0)| > Cder > 0. Let zg £ D(0,r), 
£ = 1,2,3, such that Z\ + Z2 + 23 = 0, and the minimal angle of the triangle A(zi,Z2,Zs) is bounded from 
below, and denote h = maxg \ zt\. Then, the affine map A(z) = az + /3z + 6, a, (3,5 £ C, defined uniquely by 
A(Zi) = f(ze), £ = 1, 2, 3, satisfies: 

(1) \a-f(0)\=O(h). 

(2) \p\=0(h). 

(3) D{A) = 1 + O ( ^(Jzrf^ j> with the constant inside the O-notation depending only upon the minimal 
angle of the triangle A(zi, z<i, Z2) spanned by Zi,Z2, Z3. 

(A) min fceZ |aiga i -aig/'(0) + 27r*| =G(h). 

Proof. We use the Taylor expansion (see for example, [5] page 179) of / developed around z = 0: 
(3.1) /(*)-/(0) = /'(Q)*+ Z l ' JiL]< 



27ri j c e^-z)' 

where C — <9Closure(O(0, r)) is the circle of radius r centered at the origin. We bound the reminder term 
for f(z ( ), £=1,2, 3, as follows. For z e , £=1,2, 3: 



27Ti j c e(c-ze) 



2 



h 

~ r(r — h) 



Denote, for brevity e = r ^_ h y ■ We can write three equations that characterize the affine map A(z) 
az + (3z + 5 : 

az t + m + S= /(0) + f'(0)z ( + e e ,£ = 1,2,3, 
where \e^\ < e. After rearranging: 

(a ~ /'(0)) Zi + pzj + (6 - /(0)) =e t ,£=1,2,3. 

Let us consider the matrix of the linear system: 

A = 




'We denote by Closure(Sl) the closure of the set fl. 
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Cramer's rule implies that we can bound \a — /'(0)| and \(3\ by bounding |detAf/ det A\, I — 1,2, where 
An is identical to A except that we replace its I th column with the vector (ei,£2,£3)*- A direct calculation 
shows that 

Re(z 3 - zi ) Im(z 3 - z\) 
Re(z 2 — Zi) Im(z 2 — Z\) 



det A = (z 2 - Zi)(z3 - zi) - (z 3 - zi)(z 2 - Z\) = 2i det 



Therefore, |det j4| = 4|A|, where |A| denotes the area of the triangle A spanned by zi,z 2 ,z 3 . The terms of 
the form |det ^4^| can be bounded by 

| det At\ < 6eh, £ = 1,2. 

Combining the above we get 

where in the last equality we used the fact that the minimal angle in triangle A is bounded from below and 
therefore |A| > c'h 2 . And similar bound holds for |/3|. 

Since |/'(0)| > c ( i er > 0, there exists some constant depending on Cd er such that, up-to an integer multipli- 
cation of 27r, the difference arg/'(0) — arga are of the same order as \a — /'(0)|, namely, 

min|arga - arg/'(0) + 2irk\ = 0(h). 



Finally, 



IH-|^|| \a\+0(h) 



where we used the fact that |/'(0)| > Cder > to bound a = f'(0) + 0{h) away from zero. The lemma is 
proved. □ 



Proof, (of Theorem 3.1) We construct $ 9 by sampling the uniformization map v f', that is, for Vi € V 9 , we 



set = \&(t>i), and extend <I> 9 to the whole S q by requiring linearity in each face. 

Each face fj £ F q , q > Q, is associated with some chart {f^,^}, i = ij, (used to define its conformal 
structure, see Section j2J). We denote, as above, by Vj e ,£ = 1,2,3 the (ordered) vertices of fj, and Zt = 
ifi(vj e ), £ = 1, 2, 3 their image in the local coordinates. 

We start with considering fj that are associated with interior charts, namely, charts {Qi,tpi} for vi £ 



Interior(6>). Lemma 2.3 assures existence of a constant £o > such that U p6 /.I>(<^i(p), £q) C A;. Let 
us denote by z the centroid of the triangle A {z\, z 2 , z 3 ). Let us note that the edges' length of S q are 
asymptotically of order 2~ q . In the local coordinates (i.e., Aj) the edge length goes to zero not slower than 
2~ qKi , where Ki = min{7i, 1} (i.e., depends upon the vertex angles 6i and how different they are from Oi - 
angle deficit). Let us denote h = max^ = i i 2.3 \z£ — z\ the radius of disk around z containing A (z\, z 2 , zs), and 
the conformal map ipi = ^ ° f^ 1 ■ The derivative of the conformal map, ip[, can be bounded from below over 
Closure(Ai), namely |V^(z)| _ Cder > 0, z £ A; (since (fi, and therefore ipi can be extended to a neighborhood 
of Closure (A^)). Since we have a finite number of charts and faces in F we can take Cd er to be a uniform 
lower bound for all charts. 



We now would like to use Lemma 3.2 We assume w.l.o.g. that z = 0. We set r = £07 an d note that 
T C Closurc(2?(0, 1)). Also note that all the faces in every subdivision level S q are similar to the original 
faces F of S. Lemma[2~5]therefore implies that the triangle A = A(z l7 z 2 , z 3 ) have bounded angles from below. 



Lemma 3.2 now implies that the conformal distortion of the unique affine map Aj taking A to the triangle in 

T spanned by ipi(ze),£= 1, 2, 3 is bounded by 1+0 (f^y) = l+0(2- K * q ), since h = 0(2- qK *). In addition 

this lemma indicates that |arg^(2) — arga^l = (D(2~ Kiq ), where aj = d z Aj (as-usual, the arguments are 
considered up-to addition of 27rfc). 

If we fix a domain fi C S such that Closure(f2) n V = then for fj C ft the edge length of A(zi, z 2 , Z3) is 
asymptotically 2~ q and we achieve that the conformal distortion of Aj is bounded with 1 + 0(2r q ~). Lastly, 



for this case, note that Lemma 3.2 also indicates that for sufficiently high q the affine map Aj is orientation 
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preserving and non-degenerate. Indeed since 1^(2)1 > c der > one can use Lemma 3.2 to show that for 

*,f-|&| 



sufficiently large q, detA, = \aj\ 2 — \(3j\ 2 > 0. 



Now we move to faces fj associated with boundary charts. In this case we can use Schwarz reflection principle 
to extend tpi to the union of Aj and its reflection over its straight line boundary. Note that we can bound 
the derivative of the extension away from zero also in these cases (find a new Cder > that works also for 



boundary charts). Now we are again in the situation where we can apply Lemma 3.2 (we can adjust the 
constant Eq > to also work for these charts.) 




The last case includes faces fj associated with corner charts {f2,, <fi}, i = 1,2, 3. This case can also be dealt 
with the Schwarz reflection principle, as follows. By definition of our conformal structure, the uniformization 
map in the corner coordinate charts tpi : ^ o ip^ 1 : Aj — > T, i = 1,2,3, takes the domain Aj = ^(fij) to a 
neighborhood of one of the corners of the equilateral T. Since the angle of the corner of Aj equals | which 
is also the angle of the equilateral, we can use Schwarz reflection principle to extend ipi to a conformal map 
in a neighborhood of Aj, as shown in the figure above. We first reflect w.r.t. two each of the edges touching 
the corner, and then reflect w.r.t the new line boundary (similarly to other, non-corner, boundary charts). 
The extended map is analytic (and conformal) also at the corner point because it is continuous there and 
conformal in its neighborhood. 

To finish the proof we need to show that <f> 9 are homeomorphic mappings from S to T. Note that <& 9 |/. 
is an afflne map that is composed of two affine maps. The first, denoted by Bj, takes fj in S q to the 
coordinate chart's triangle A = A(z\, z 2 , z 3 ), and the second, Aj, maps A to A^^-j), ^(vj 2 ), ty(vj 3 )). The 
first affine map is orientation preserving and non-degenerate (and with bounded conformal distortion) for 



sufficiently large q due to Lemma 2.5 The second affine map Aj was proven above (by showing that the 
determinant of its Jacobian is positive) to be orientation preserving and non-degenerate as well for all faces 
and for sufficiently large q. Therefore, <i> q jj. is orientation preserving non-degenerate as a composition of 
two orientation-preserving and non-degenerate affine maps. Since <I> 9 also maps the boundary of S bijectively 
onto the boundary of T (by definition of J- s ")), it is a global bijection (see [T7] for a detailed proof). Since 
it is simplicial, its inverse is also continuous, a fact which makes it a homeomorphism. □ 

A comments is in order. One can build the charts based on an already subdivided version of the surface S q , 
q > 0, and achieve better estimates of the conformal distortion of the maps $ q . For example, if one takes 
the charts defined by S 1 , the vertices Vi € V 1 \ V are all flat, i.e. 0i = 2ir, and therefore triangles associated 
with the corresponding charts will have conformal distortion of order 1 + (D(2~ q ). 
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4. Approximation 

Our next step is showing that if we are able to put our hands on a series of simplicial maps with discrete con- 
formal distortion that goes uniformly to one, then this series has to uniformly converge to the uniformization 
map Note that this result applies to any such series, not only the ones built in Theorem |3.1| 



Theorem 4.1. Let \& : S — > T be the uniformization map of the polyhedral surface S. Let {$ 9 } g >Q C J~f?, 
for some series K q = 1 + 0(2~ qK ), where n > is some constant. 

Then, $ 9 o if -1 : 7~ — > 7~ converges uniformly in any compact subset of Interior(T) to the identity map 
I, : 7 > 7 . 

Proof. Define g q = $ 9 o ty' 1 : T -> 7~. 

We start by showing that g q is TT-quasiconformal with some global distortion bound K > 1. 
First, since \& is an orientation-preserving homeomorphism, "P -1 is orientation-preserving homeomorphism. 
By assumption $ 9 is an orientation-preserving homeomorphism and so g q is an orientation-preserving home- 
omorphism. Let us observe that g? = g q |*(i n tcrior(/j)) can be written as 

g q = $« o ip^o ipj k(Interior(/j))i 
(*) 

where i — ij is the chart associated with the face fj £ F q . The part that is marked with (*) is a conformal 
map over \& (Interior(/j)) (by the definition of the uniformization map). Next, recall that tpi consists of a 
composition of a rigid maps and power map h{z) = z li , which is conformal except at the origin. Therefore, 
ip^ 1 is conformal over the domain ipi o vp^ 1 (vp (Interior(/j))) = (pi (Interior (/j)). Lastly, the <fr 9 part is 
defined over the set ip^ 1 o tpi o vf^ 1 (v? (Interior(/j))) = Interior(/j) and can be written as a composition 
of two affine maps (as before) the first, denoted by Bj, takes fj in S q to the coordinate chart triangle 
A = A (ip^Vji), (fii(vj 2 ), <fi(vj 3 )), where Vj t , I = 1, 2, 3 are the (ordered) vertices of face The second affine 
map, Aj, maps A to A(^E , (wj 1 ), ^(«j 2 ), ^(wj 3 )). Bj is i^'-quasiconformal, for some K' > 1 independent of q 
by Lemma 2.5 is _ftT"-quasiconformal, again independently of g, due to the assumption that $ 9 G ^ r ^ 9 , 



= 1 + 0{2- qK ). Then, g] is K = ^'^"-quasiconformal, for all faces fj £ F q and all q. Denote the 
domain 

D = U/ j£F «* (Interior (/,)) C Interior(T). 

The set T\D, which consists of the image under ^ of the edges E 9 of S q is a union of analytic arcs. By 
extension theorem for quasiconformal mappings (see e.g., Theorem 8.2, page 42, and Theorem 8.3, page 45, 
in [TS]) the conformal distortion of g q over the whole T is K again. Since bounded if-quasiconformal maps 
form a normal family there exists a subsequence g qk that converges locally uniformly to a if-quasiconformal 
map or a constant (see Theorem 3.1.3, page 49 in [5]). Denote the limit map by g. Since g q fixes the corners 
of 7~, g is i^-quasiconformal. Next we show that g is a conformal map. It is enough to prove that g is 
1-quasiconformal (see e.g., Theorem 5.1, page 28 in [15]). For that end, take an arbitrary domain O C T 
such that Closure(O) n ^ (V) = 0. Now, we look at g q *"° = g qk |on*(intcrior(/ J ))- Similar to above, it consists 
of a composition of conformal maps (ipi o vp^ 1 and the inverse of the power map h(z) = z 7i ), and an affine 
map that can be represented as composition of two affine maps Bj and Aj. The conformal distortion of 
Aj equals T)(Aj) = 1 + (D(2~ qK ) by assumption. The affine map Bj is defined by sampling the power map 
h(z) — z li , i = ij at the corners of fj (the origin is placed at vertex vi £ V) and extending linearly. Since, 
by the intersection with the domain O, we sample h(z) at some bounded (from below) distance from the 



origin (vertex Uj), Lemma 3.2 indicates that D(£>j) = 1 + 0(2 q ). This implies that the conformal distortion 



of gj k ' is 1 + 0(2 qK ). As before, extension theorems for quasiconformal mappings implies that g qk \o is 
I + 0(2 -9K )-quasiconformal. Taking arbitrary large Q and observing the series tail 4> 9 for q > Q implies 
that g \o is 1-quasiconformal (for relevant convergence theorem, see for example Theorem 5.2 in [15] page 
29). So we get that for arbitrary domain OcT such that Closure(O) n ^ (V) = 0, g \o is conformal. This 
means that g is conformal in Interior(T) \ ^(V). However ^(V) is a set of discrete points and therefore 
standard extension theorems of conformal mappings imply g is conformal in the whole Intcrior(T). Since g 
fixes the corners of the equilateral it has to be the identity map since it is the only conformal bijection from 
T —> T that fixes the three corners. 
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Since every infinite subsequence of {g q } has a subsequence converging locally uniformly to the identity, g q 
has to converge to the identity map locally uniformly. 

□ 



5. Convex Spaces of Simplicial Quasiconformal mappings 



In this section we develop the convex spaces of simplicial quasiconformal mappings, prove the main theorem of 
this paper (Theorem and develop two algorithms for approximating the uniformization map \P : S — > T. 

Let us recap what we know so far: 

(1) Uniformization: there exists a unique conformal homeomorphism of the polyhedral surface to the 
equilateral : S — > T taking three prescribed boundary vertices Wi,U2,i>3 to the corners £1,^2,^3 of 
the equilateral T ■ 

(2) Feasibility: for sufficiently large q, Ff^^-i*) 1S n0 ^ em Pty> where k — min{Kj}, = min{l,7i}. 
Moreover, J 7 ^ 9 is not empty where K = (K\, Ki-p«\), and Kj — l+0(2~ qKi ) for fj e F q associated 
with the chart i = ij. 

(3) Approximation: any series {$q} g >Q, &q € •^x-\-Q(2-i^y wnere k > is a constant, converge to the 
uniformization map $ on compact subsets of Interior(T). 

These observations motivate the meta-algorithm shown in Algorithm [T] for approximating the uniformization 
map 

Algorithm 1: Approximation of the uniformization map. 

Input: Polyhedral surface S = (V, E,F), with three positively-oriented boundary points v\, V2,va G V 

Subdivision level q > 2 > 

Constant < c < 1 
Output: Quasiconformal simplicial map $ 9 

// Set the charts 
forall Vi £ V do 

_ Calculate 9i , 7$ , Ki ; 

// Subdivide 

Subdivide S q times to get S q = (V q ,E q ,F q ); 

II compute discrete conformal structure 
forall fj eFdo 

Associate fj with a chart i = ij] 

Choose a coordinate system on fj such that Vi is placed at the origin; 
Map the vertices Vj e ,£ = 1, 2, 3 of fj by z 1— > z 7i , and set as Aj; 
_ Kj <- l + 2- cqK >; 

K <— (K 1 ,K 2 , ..,if|F9|); 

// Approximate 
Find a map $ 9 g ; 

Return 



The correctness of this algorithm can be ex plained as follows: since < c < 1, the series 2~ cqKi i approaches 
zero slower than 0(2~ 9Ki j ), Theo rem 3.1 guarantees that for sufficiently large q the set with Kj — 
1 + 2~ cqKi j is not empty. Theorem 
uniformly. 



4.1 



then implies that $ 9 will converge to the uniformization map locally 



At this point we are left with the problem of finding an element in the non-empty set J-^ . This is a non- 
trivial task since J-^ is a non-convex set and hence finding an element in this space is a non-convex problem. 
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Nevertheless, as we explain next, we can identify a collection of convex subsets of T^ 1 the union of which 
spans T^ . 

The space T^ contains orientation-preserving homeomorphic simplicial maps $ 9 such that D d ($« \ f 
) = T)(Aj) < Kj, where Aj is defined by mapping the triangle Aj = Aj {z\, Z2, 23), zi = ipi(vj e ), 
£ = 1,2,3 (i.e., the triangle formed by the vertices of fj written in local coordinates) to the triangle 
A($ 9 (d :( - 1 ), <fr q (vj 2 ), & q (vj 3 )) (i.e., the image of fj under Next, we write down characterizing equations 
lor ./•£'. 

We start by examining a single face fj g F q . We will work with its image in the local coordinates as given 
by its associated chart {Qi,(pi}, i = ij. An affine map A — Aj mapping A-,- to the complex plane can be 
written (using standard complex notation) as 

A{z) = az + j3z + S, a,p,5eC. 

In what follows we will only consider the linear part of A, namely az + /3z. We will abuse notation and denote 
A(z) = az + /3z, the translation does not affect conformal distortion nor the orientation and therefore it will 
be safe to ignore it for now. We ask that A is orientation-preserving and D(„4) < K, for some arbitrary but 
fixed K > 1. The conformal distortion of A can be calculated using the formula (see e.g., [T], page 11) 

(5-1) V(A) = g±j§ < K. 

\a\ - \/3\ 

And the Jacobian of A, 

(5.2) 3(A) = \af-\f3\ 2 . 

A is an orientation preserving homeomorphism if and only if J(*4) > 0. Therefore, necessary and sufficient 
conditions for the affine map A to be if-quasiconformal are 

(5.3) D(^) < K 

(5.4) 3(A) > 



We will write the above conditions in a more convenient form. First, using (5.2), (5.4 1 is equivalent to 
(5.5) |a| > \p\. 



Second, using (5.1), (5.3) can be rearranged, taking into account that (5.5) implies a/0we get 



(5-6) o^O 

Let us denote the collection of orientation preserving planar linear transformations with conformal distortion 
K by Tk- The collection Tk can now be thought of as a subset of C x C where a pair (a, 13) G C x C 
represents an element (without translation) A £ J-r, namely, A(z) = az + (3z. 

In this parameterization we characterize the maximal convex subspaces of Tk- As we prove next there is 
one parameter family of maximal convex subspaces Tk.t C Tk C C x C defined by: 



(5.7) |0| < |^Re(V iT a), a^O 

From this definition we see that J~k,t ar e convex subsets of Tk (this can verified directly by taking convex 
combinations of elements in Tk,t)- These convex spaces span Tk, 

Tk = [J Tk, t - 

tG[0,2tt) 

Lastly, Tk,t are maximal convex subsets of Tk- We say that U C Tk is a maximal convex subset if every 
convex set V C Tk that contains it U C V has to be equal to it, U = V. 

Lemma 5.1. Tk,t, t S [0,2tt) are maximal convex subsets of Tk- 
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Proof. Let U C Tk be a convex subset such that Tk,t £ U. Let (a* , f3*) £ U\ Tr,t- By definition (5.7) 
Tk.t contains the linear map (e 1T £^\/3*\, f3*). Now since U is convex it contains the convex combinations 

(a(A),/3(A)) := ((1 - A) e lT f^|/3*| + Aa*,/?*) 6 17. 



Since (a* , (3*) £ T K we have that 
Using this inequality, 



KA)| 2 < (|^) 2 |/31 2 (1-A 2 ) + K|A 2 



^ + n 2 ,,„2 , , 2 f. „ /^ + i\ 2 |2 N 



Hence, for sufficiently small A we get that 

K + 1 A' + f 

KA)|<^|/r| = -^|/3(A)|, 

which is a contradiction with the fact that [7 C Tk- O 

The subset Jjf, T C Tk contains affine maps A(z) = az + f3z + S £ J 7 ^ with distortion D(_4) — k < K that 
in addition satisfy 

/ x -i fK + lk-V 

(5.8) |arga-r|<cos ^^y^y 

As usual, the argument in the l.h.s. should be understood up-to addition of k2?r, k E Z. To see this, let 



A{z) — az + j3z + 5 have conformal distortion D(_4) = k < K, and assume it satisfies equation (5.8) then 
we need to show that A G J~k- Indeed, we have 



> 



Re(e 1T a) = \a\ cos (arg(a) — r) 
K + lk-1 
K-lk + 1 

- 

where in the last equality we used the fact that D(A) = k. These inequalities also show that the angle in the 



r.h.s. of (5.8) cannot be enlarged while A £ Trt- This shows that (5.8) is necessary and sufficient condition 



for affine map A with distortion D(^4) = k to belong to Tk.t We proved 

Lemma 5.2. The convex space J~k,t contains all affine maps A with conformal distortion k — D(A) < K 
that their rotation (i.e., &rgd z A) satisfies equation (5 



A corollary of this observation is that there does not exist a convex subset U C J~k such that it contains 



all affine maps A with conformal distortion k — D(_4) < K with larger angle limit than described in (5.8) 



Indeed, if there was such a convex space it would be a superset of J-k,t which will contradict the maximality 



of J~x T shown in Lemma 5.1 We summarize: 



Corollary 5.3. There is no convex subset of the K -quasiconformal simplicial mapping space, U C Tk, that 
contains all the affine maps A with conformal distortion D(A) < K and larger rotation angle limit than 
J~k.t for any r £ [0, 2ir) . 

Intuitively, the maximal convex space Tk,t restricts the rotation angle of its member affine transformations 
A (i.e., arga) to be the maximal possible around a prescribed rotation by r radians. For perfect similarities 
A(z) — az the rotation angle range is (t — f ,t + ^), while for maps with conformal distortion k = K only 
rotation by exactly r is allowed. Affine maps with intermediate distortion values 1 < k < K will be allowed 



rotations in between these two extremal cases, as expressed in equation (5. 
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We now move to the general triangulations S q , and their corresponding spaces J-^ 1 of homeomorphic simpli- 
cial maps : S q ->■ T where D d ($<? \ f .) < Kj, K = (K u .., K\ F «\). Our goal is to find an element in Tff . 
We start by formulating a set of equations that characterize exactly. We derive necessary conditions 
for $ 9 G ) an d later show that they are also sufficient. Let <& g G ' then, over each face fj S F', $ j^. 
is an affinc map, in the local coordinates it has the form 

Aj (z) — a.jZ + PjZ + Sj . 

Also set Ui = G C. By definition we have D(A,) = Kj < K. That is, Aj should satisfy the equation 

Ki-1, 



(5.9) 



1 



Furthermore, since $ 9 is continuous, each affine map Aj should map the vertices of its face fj to U{. Namely, 
if we denote by ji, £ = 1, 2, 3, the indices of the vertices of face fj, and Zj f = ifi(vj e ) their complex coordinates, 
then 



(5.10) Aj {zj t ) = a jZjt + (3 j z je + Sj = u jt , 

Note that these are homogeneous linear equations. 



1,2,3. 



Another set of necessary conditions is related to the boundary conditions. We want the boundary dS q to 
be mapped to the boundary of the equilateral dl~. This can be achieved by first forcing the corner vertices 
w ij v 2, to be mapped to T's corners ^1,^2,^3- For each face fj that contains one of the corner vertices vi, 
£ = 1, 2, 3 we add an equation of the form 



(5.11) 



Aj(zi) = Q.jZt + fa z e + S 3 = t e . 



Second, the rest of the boundary vertices Vi G dS q are mapped to the edges of T . This can be achieved 
by ordering the vertices that are supposed to be mapped to one edge of T, say {t\,t^), w.r.t. the positive 



orientation: Uj 
equations 

(5.12) 
(5.13) 



v 1 ,v 



121 '") u l>n — 1 ' l n 



v 2 - Now introduce new real variables A2,...,A„_i and set the 



< A 2 < A 3 < ... < A n _i < 1 
He = (1 - A^)ti + Af* 2 , £ = 2,..,n—l. 



We can use < instead of < in equation (5.12) since (5.9) implies strict inequality automatically. 



To recap, we derived necessary conditions for $ 9 to belong to J 7 ^ . Le t us sh o w tha t they are also sufficient, 
that is, every piecewise affine map that satisfies equations (5. 9), (5. 10), (5. 11), (5. 12), (5. 13) belongs to J-^ ■ 



Proposition 5.4. Let<fr q be a piecewise affine map overS q that satisfies equations (5.9), (5. 10), (5. 11), (5. 12), (5. 13) 
Then, $9 e jrS 9 . 



Proof. From equations (5.10) we get that <i> 9 is continuous. Equations (5.9) imply that over each face fj, $ is 



an orientation-preserving affine map with conformal distortion bounded by Kj . This implies that $ : S q — > T 
is a bijection (see [17j for a proof). Since <& 9 is an orientation preserving piecewise-affine continuous bijection 
it is an orientation preserving homeomorphism. Therefore <& q G ■ Q 



lj w e can find an element in J-^ by simply looking for variables {ui} , {aj, fij, Sj} ,{Xi} 



Going back to Algorithm 

that satisfies equations (5.9|),(|5.10|),(|5.11|),(|5.12|,(|5.13|). In fact, equation (|5.10|) can be used to express the 



variables {a 3 ■ , @j, Sj} as linear combinations of the variables {u^ and therefore eliminate these variables from 



the equations. All the mentioned equations are convex equations except equations (5.9). Since the bounded 



conformal distortion equations (5.9) are not convex, finding a feasible solution <& 9 is not a convex problem. 



However, we can replace equation (5.9) with the convex conditions (5.7) (describing the maximal convex 
subset) and achieve a convex problem. That is, for each face fj G F 9 , we need to assign an angle Tj G [0, 2ir), 
and replace (5.9) with the equation 



(5.14) 



< 



K1 



a, f 0. 
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Let us denote by J-^.t > T — ( r i> t 2j — j t |F9|) G [0, 27r)l F9 l the convex space denned by equations ( |5.14 1,(5.10 1, 



(5.11 1,(5.12 1,(5.13 ). By construction and Proposition 5.4 T§^ T C . In Algorithm [l] we set K 



1 + 2 _C9Ki 3 , where ij is the index of the chart associated with face fj. We proved in Theore m |3. 1| that there 
exists a simplicial map, let us denote it <£> 9 , such that the discrete conformal distortion D d ($ 9 that is 
D(Aj) equals 1 + 0{2~ qKi ' '). Denote AAz) — oijZ + (3jZ + S, and Tj — argSj. Equation (5.8) now implies 
that J'kV wm contain <J> 9 if 

■ 2 + 2 -cq Klj QQ-IKH) 

for all fj G F 9 , and the r.h.s. converge to ^ as <? — > oo (remember that the constant inside the 0-notation is 
set per chart and there is a finite number of charts). Therefore, for any e > 0, there exists sufficiently large 
q such that $ 9 G J~^ T , as l° n § as \ T j ~%\ < § — e - ( as before, the arguments should be understood up-to 
addition of k2ir, k G Z) We have proved: 



\ta — tA < cos 



3.1 



Lemma 5.5. Lei <&j : S q ^ T be the simplicial map, the existence of which is set by Theorem 
its restriction to face fj in local chart coordinates by Aj(z) = cijZ + fljZ + S, and set Tj = argSj. 
For arbitrary e > there exists sufficiently large q such that if T = (r 1; - . , 7"i jre i ) satisfy 



Denote 



max mm r. 



fc27r| < 



then, 



& G Jf> 



We are now in a position to prove the main theorem of this paper: 



Theorem Let ^ : S — > T be the uniformization map of a disk-type polyhedral surface to the equilateral 
T taking three prescribed boundary vertices of S to the corners of f . Let S q be the q th -level subdivided 
version of S. 

If, for an arbitrary but fixed e > 0, the argument of 4 r ' is known up to an error o/±(^ — e), then one can 
construct a series of non-empty convex spaces U q = {$> q } C J- sq of simplicial maps of S q such that: 

(1) Every map <& q G U q is K-quasiconformal (QC) homeomorphism that maps S onto T, with some 
constant K independent of q. 

(2) Every series {& q } q >Q 7 where $ 9 G U q , converges locally uniformly to the uniformization map ty. 
That is, <I> 9 o vfi^ 1 : T — > T converges uniformly in any compact subset of Interior{T) to the identity 
map Id ■ T — > T '. 



Proof. Let e > be arb itrary but fixed. Let $ 9 : S q — > T be the simplicial map, the existence of which is 



assured by Theorem 
1, 2, 3, and A 



3.1 



Denote, as before, Aj{z) = o.jZ + /3jZ + 5, and Tj — argSj-, zg = ipi(vj e ), i — ij, I 



A(zi, Z2, Z3) the face fj in the local coordinate chart, and the centroid z = | {z\ + z^ + Z3). 
Let (by Theorem 3.1 again) Q > be sufficiently large such that for all q > Q, fj G F 9 , 



(5.15) 

Now we set U q 
T = (n,--,T|p5i| 



mm 

fcez 



Tj - arg 



oc Pi 



(?) 



fc27T 



< 



|F«M 



Kj = 



1 + 2 C9Ki , for some constant 1 > c > 0, 



where K = \K\,..,K\ 
is defined by equations (|5.14|),(|5.10|, (|5.11|),(|5.12|),(|5.13|), and Proposition 



5.4 



implies that -FjlV C Ff? , in particular, it contains only X-quasiconformal homcomorphic simplicial maps 
from S — > T , for some sufficiently large but constant K. By assumption we know the argument of \&' up-to 



| — s and so we take some Tj such that 



(5.16) 



arg 



k2ir 



< 
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By triangle inequality, (5.15) and (5.161 imply 



mm \Tj 

fcez J 



k2-rt\ < 



Now we want to use Lemma |5.5| Set Q > sufficiently large as required by Lemma 5.5 (we could have set 



Q > sufficiently large a-priori) and we get <& g £ -^kt' ano - so ^ ^ s n °t em Pty> f° r all g > Q. Now take an 
arbitrary series {<f> q } q>Q , $ 9 £ T^ T , where Kj = 1 + 0(2" c « Ki ) = 1 + 0(2-1"), k = cmm{/ti} > 0. The 



convergence of the series $ 9 to the uniformization map is implied by Theorem 4.1 



□ 



An immediate consequence of Theorem |l.l| is that if we have approximations Tj of the argument of Vl/' up- to an 
error of ±^ then we can simply solve the convex feasibility problem <£> 9 £ -^kV t° achieve an approximation 
of the uniformization map. In particular, this can be achieved by solving the following convex minimization 
problem: 



(5.17) 



mm e 
s.t. 



and 



< 



A, 



A, 



1 



-Re 



eqs. (5.10), (5.11), (5.12), (5.131. 



Next we suggest two algorithms on how to set r without any prior knowledge: the first is exhaustive, guaran- 



teed to find the feasible r but requires solving (5.17) very large number of times (exponential in the number 



of faces of S). This algorithm is guaranteed to work, but its computational complexity makes it impractical 
for applications. The second algorithm is an iterative one that start from arbitrary r° and construct a scries 
of t 1 ,!- 2 , ... until finds a feasible space J-^ T - The second algorithm is much more efficient than the first 
one, however lacking a proof that a feasible space is always found (for sufficiently high q). Nevertheless, once 
found a solution, all the theoretical guarantees apply and the convergence and approximation results apply. 
In practice this algorithm works well as we demonstrate at the end of the paper. 



5.1. First Algorithm (exhaustive). We can exhaustively search for a feasible r by testing for each Tj 
three angles, say Tj £ {0, X'Tf} - This way we are guaranteed to find an assignment r = (ti, ..,t|f«|) 
satisfying the condition of Theorem |1.1| However, testing all such assignments requires solving (5.17) for 
3' F '! convex spaces J~^ T , which is computationally daunting. However, it is clear that the argument of the 
derivative of the uniformization map cannot change arbitrarily in small areas, and indeed, as expected, we 
will show that the number of convex spaces needed to be searched is at-most of order 3' F ', which is still high 
but at-least independent of the subdivision level q. 

It is possible to prove this by developing an approximation argument for conformal maps, using the ap- 



proximation Lemma 3.2 (observing how the argument of the derivative of the affine maps approximate the 
arguments of the derivative of the conformal map). However, we will pick a different route which we believe 
also gives an interesting intuition regarding quasiconformal simplicial maps. Intuitively, the claim is that the 
argument change in simplicial quasiconformal mappings is linearly (and simply) bounded by the conformal 

distortion bound. Denote by dgi ^e 1T , e 1T J = min^gN |r — r' + 27rfc| the distance on the circle S 1 , then, 

Lemma 5.6. Let fj,fj> be two planar faces sharing an edge. Let Aj,Aj> be two planar affine maps, coin- 
ciding on the common edge and have K conformal distortion bound, that is D(Aj) < K, and D(Aj>) < A. 
Then 



' s 1 



\d z A 3 y \d z A 



< sin 



> K ~ 1 
'K + 1 



Proof. Let e& denote the common edge of the faces fj,fj>, and let = (UjUj/), Vi,Vi> its two vertices. By 
continuity 

Aj{v v - Vi) = Aj>(vi< - Vi). 
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Figure 3. Illustration for the proof of Lemma 5.6 



Plugging the expressions Aj{z) = ajz + j3jZ + Sj, Aj(z) = ayz + /3j/z + Sjr for Aj,Aji we get 



aj(vi> - Vi) + /?.,•(«»/ - Vi) = a 3 i{v v - v t ) + /3,<(tv - v t ), 



and so 



Taking the absolute value and dividing by |iv — Vi\ we get 

\<*j-<*j'\ = I ft' 

SinceD(ylj) < if, we have that |/3j | < | ctj | , and similarly for Ay. Let us assume w.l.o.g. that \oij \ > |<x,-< 
Therefore, 

if — 1 ,, ls K- 1, , 

k + k' < 2 — — rla,-' 



if + 1 



Dividing both sides by |ay|, we get 



1 



< 2- 



( if - 1 



It is not hard to check now (see Figure [3]) that 



arg ■ 



< sin" 1 | 2 



if- 1 
if + 1 



□ 



Lemma 



5.6 



indicates that for $ 9 e J 7 ^ > ^ 



to the same chart {fii, tpi}, 



1 + 0(2 qKi ), a. pair of adjacent faces fj, fj> C f2j associated 



O (2 



where we used a Taylor expansion of sin _1 (-) around zero. This is still not enough since to traverse from one 
face to another in the q th subdivision level inside fii we will need 0(2 q ) edge crossings, and hence we only 
get that the difference between different arg o.j for fj C % is only bounded by 0(2 q{1 -^) and since Ki can 
be strictly smaller than one, this is still not a constant. We therefore need to refine our argument. We will 
show that the argument change in a chart Q$ can be bounded independently of the subdivision level q. 

To do that we note that in Theorem El] the estimate D d ($« \ f .) = D(A,-) = 1 + 0(2^ qKi ) (remember that 
Ki = min {1, 7i}) can be refined by replacing it with 

D d (<F \ fj ) = n(A j ) = l + 0(h j ), 

where hj = max^i^^ \zg — z\, z = \ {z\ + Z2 + Z3) is the centroid of the triangle A(zi, z%, 23), and zg = 
ipi(vj e ) are the images in the local coordinates of the vertices of the face fj, and the constant in the de- 
notation is independent of which face fj we choose in chart {Q,i,tpi\. Now, take two faces fj,fj> C Hi in 
chart i = ij. W.l.o.g. we can assume they both belong to one / 6 F, that is fj, fj> C /, otherwise we can 
use the following argument few times but not more than half the valence of the vertex Vf (defining the chart 
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i). We can also assume that fj is the face touching the vertex Vi, otherwise we can bound the argument 
change by the sum of bounds when traversing from fj to the corner face (touching Vi) and from the corner 
face to fj>. 

Starting from (the corner) fj we traverse to any fj/ in the manner depicted 
in the inset figure (red triangles represent the path connecting the face fj, in 
yellow, and the face fj/, in green). Denote the faces we traverse by fj — > fj 1 — > 
f j2 fji. We would now like to bound 0(hj) + 0{hj t ) + 0{h n ) + 

... + 0(hj>). Since, as mentioned above, the constant in the O-notation is 
independent of the choice of face in the chart we actually need to bound hj + 
hj 1 + hj 2 + ... + hj' . This will be done by using the binomial formula. 

The first thing we do is to bound h for a general face fj" C /. We will then use 
this bound to bound the sum of h's for our path. To bound h = rnaxi \ zi — z\ it 
is enough to bound the length of one of the edges of the triangle A(zi,Z2, Z3), 
where zg = ipi(vj"), and Vj",£ = 1,2,3 are the vertices of the face fj". Let us 
denote by m the generation of fj" counted from the corner, as shown in the 
inset . 

Denote by Wi,W2,W3 the corners of fj". A(wi, W2, W3) is similar to /, and 
w.l.o.g. we assume Wi is the vertex corresponding to /'s corner under the 




m = Q m = \ m = 2 ... m = 2 q 



similiarity (see inset figure for illustration), 
edges of A(zi,z 2 ,z 3 ) is 



w4 



Then the length of one of the 



+ ('«i + (w 3 - wi)) 7l | 



(5.18) 



where we set \x = 



W3—W1 



= Krii-a+^i, 

Now using the binomial expansion 



Wi 



U'l 




(l + /ip = l+7i// + 



7t 
2 



3 



= 1 + jifi + i?, 



and the reminder term R can be bounded by noting that w.l.o.g. 7,; < 1 (otherwise the edge length is 0(2 q ) 



and there is nothing to prove) and therefore by induction 



li 
k 



< 1, for k > 2, then, 



1*1 < H 



Plugging back in (5.18) we get 



,,Ti I - 



1 - h 



|1-(1 



< (iM + 



7iM + *)l 

IH 2 



1 



I Ml 



It is not hard to see that there exists constants c\,C2 > such that \w\\ < C\m2 q , and < C2m , and 
therefore 



(5.19) 



1 1 < c 2' q ^ 



where c > is some constant. Now we use (5.19) to bound the sum of edge lengths in our path described 
above, 



hj + h n 



< 



c2 -1'H ^ w 
m— 



74- 



< c2- q ^ (C»(2 TO ) + 2 n ) 

= om, 
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where we used the fact that m < 2 q , and that V 



This means that we can consider only r = m, .., t]f«i) that are defined on a subdivision level of, say, 
q' = [log 4 <7]. That is, for subdivision levels q = 2,3,4,5,6,..., we take the assignment vector r to be 



piecewise constant on faces of S q , q' = 1,1,1,2,2,... (respectively). Since 



O (4« 
If 9 ' 



= O(g), the 

number of assignments r we need to check to find a feasible map is of order O I 3l* " M = O (3«). The 
arguments abobe imply that for sufficiently large q there exists some assignment constant over faces of level 



q' that will lead to a feasible convex problem (5.17) 



A comment is in order: as the argument of conformal maps is harmonic, it is very likely that the argument 
change will be slow even between neighboring charts in S (not even the subdivided versions). Therefore, it is 
possible to reduce the complexity of the algorithm for many polyhedral surfaces S by defining assignments 
by spreading small number of point seeds over S, dividing S into sub-areas of "constant" argument and 
transporting the argument using the chart. Obviously, the way the arguments will be transported will change 
the transported values, but nevertheless this type of argument will result in an algorithm with complexity 
connected only to the geometry of S and its uniformization map to T (and not to its flat "tessellation"). 



5.2. Second (iterative) algorithm. Here we suggest a practical and simple algorithm for finding $ 9 € 
J-jl . In a nutshell, this algorithm iteratively adjusts r based on the result of the previous iteration. Although 
we do not have a proof that this algorithm will always find a feasible solution (we will leave this issue to 
future work), it seems to do well in practice, and in the case it does find a feasible solution it is guaranteed 
to be an approximation of the uniformization map, and has all the theoretical properties of the exhaustive 
algorithm. 

The idea of the algorithm is very simple: start with some arbitrary angle assignment t = (rf , T|p g |) and 
solve the convex feasibility problem <I> 9 € -^kt ' wnere K = {K\, ..,Kipqi), and K 3 ■ = 1 + 2~ cgKi . This is 



done by solving the cone programming problem (5.17|. If e < a feasible solution was found and we are 



done. Otherwise set rj = arga^, where ctj is, as before, d z Aj. The outline of the algorithm is provided in 
Algorithm [2] 

We have implemented this algorithm in MATLAB and used it to approximate uniformization of polyhedral 
surfaces. One useful variation of Algorithm 2 is that once a feasible solution is found € -^kV- ^ ^ s 
possible to further improve the result by pushing the conformal distortion per face further down using the 
following, now feasible, problem: 



(5.20) 



min 




s.t. 






< 




< 


and 





Ki + 



|Rc(e i7 >ay) + Ej , 



eqs. (5.10), (5.11), (5.12), (5.13) 



These equations take the role of eqs. (5.17). We now iterating with these equations, keep updating r and 
resolving until J2j £ j is n °t decreasing. 



Another comment is that in implementing Algorithm [2] it is possible to replace the cone conditions in (5.17) 
(and similarly in (5.20)) with linear constraints to end up with standard linear programming problem: 
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Algorithm 2: Approximation of the uniformization map (iterative algorithm). 

Input: Polyhedral surface S = (V, E, F), with three distinct boundary points Ui,i>2)t>3 € V 

Subdivision level q > 2 > 

Constants < c < 1 
Output: Quasiconformal simplicial map $ 9 

// Set the charts 
forall Vi £ V do 

Calculate Oi,~fi,Ki; 

II Subdivide 

Subdivide S q times to get S q = (V 9 , E 9 , F 9 ); 

// compute discrete conformal structure 

forall fj e f do 

Associate fj with a chart i = if, 

Choose a coordinate system on fj such that Vi is placed at the origin; 
Map the vertices Vj e , I = 1, 2, 3 of fj by z H > z 7i , and set as Aj; 

_ Kj 4r- l + 2- c « K - 

K^{K X ,K 2 ,..,K^); 

II In the following: find a map <f> 9 e jf' 
Initiate r arbitrarily: € [0, 27r); 
while e is sti// decreasing do 

Solve ( [5171 ); 

Update r: Tj- ■<— arga^ ; 

if e < then 

L Return 

else 

Return "Did not hnd a feasible element in J-^ ■ "; 



(5.21) 



s.t. 



|oo 



and 



< 



: Re 



1 



(e-^aj) 



1 Kj 

eqs. (OO), (pull), (pU2j), (I5~l3|, 



where we denote |j9|oo = max{|Re/3| , |Im/?|}. It is shown in [T7] (see Proposition 4.2), that this formulation 
will be equivalent to the cone condition as Kj —> 1. 

Figures [T] and [4] show two examples of approximations of the uniformization maps calculated with our iterative 
algorithm. The first row shows the subdivided series of meshes S° = S -< S 1 -< S 2 -< 5 3 . The second (and 
third row in Figure [I]) shows a checkerboard texture mapped to the polyhedral surface by the inverse of the 
simplicial approximations <1> 9 . The bottom row shows the homeomorphic image under <i> 9 of the triangulation 
S q onto the equilateral. Red color indicates conformal distortion. Note how the approximation improve as 
the refinement level increases. 
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s s 1 s 2 s 3 

Figure 4. Approximation of the uniformization map of the Igea polyhedral surface. 

Appendix A. 



Lemma 2.5 Let T = A(£,0,?7) be an isosceles triangle (\£\ = \rj\) and denote the angle 9 — Z(^,0,?y). 
Further let T q be the q th level of regular 1-4 subdivision of T . Denote by h(z) = z 1 the power map, and 
assume that \j~\0 < ^, and that 9 < 60.4°. Then, the simplicial maps h q defined by sampling h(z) over 
the vertices of T q and extending by linearity are homeomorphisms that satisfy D(h q ) < K for some K > 1 
independent of q.. 



22 



YARON LIPMAN 



Proof. We will denote by x q '■ T — > C the simplicial map that is defined by sampling h{z) = z 1 at the vertices 
of T q and extending linearly. In particular, x q w iU map affinely the triangle A = A(zi, Z2, 23) C T q to the 
triangle A(z 7 , z% , z 7 ). We need to prove that x q are &H quasiconformal maps with a universal bound on 
their conformal distortion (a bound independent of q). 

We first prove that it is enough to show that for sufficiently large Q > 0, x q q > Q is a homeomorphism. 
Then we show that such a Q indeed exists. 

Let us denote the set Tf ar C T to include all the points z £ T such that if we write z = Ai£ + A2f? then 



3.2 



indicates that there exists some 



Ai + A2 > \. We also denote T near = Closure(T \ Tf ar ). Lemma 
finite level of subdivision Q, such that for all subdivision levels q > Q, the map \ q restricted to the set of 
triangles in T q contained in Tf ar is quasiconformal with bounded conformal distortion of, say, K = 2. That 
is, D(x 9 \r far ) <1 ioi q> Q. Now our assumption is that for sufficiently large Q, x q is homoemorphic for 
q > Q. In particular we can take Q > Q. 

Now, x is homeomorphic, and therefore has to be quasiconformal since it maps finite simplicial complex 
(build out of finitely many triangles). Denote D(y^) = Kq. Since h(az) = (az)' y = a' y h(z), and multiplying 
by a complex number keeps distortion and orientation unchanged we see that the conformal distortion of 
•^-Q+i over triangles in T® +1 contained in T near is also bounded by Kq. Therefore, the conformal distortion 
D(x Q+1 \r naar ) < K Q , while we already know that D(x Q+1 |r /ar ) < 2. Therefore 

D(X Q+1 |t„„J <max{^ Q ,2}. 
Continuing with induction we prove that 

D( X 9 |t„„J <max{X Q ,2}, 

for all q > Q. 

We are left with the task of proving that x 9 , 9 > Q are all homeomor- 
phisms for sufficiently large Q. We will refer to the inset image depicting 
an example of A(^,0,ry) for illustration. It is enough to show that ev- 
ery triangle in every subdivision level T q is not flipped by x q ■ Indeed 
since the boundary polygonal of T q is mapped to another simple poly- 
gon, the fact that all inner triangles maintain their orientation implies 
homeomorphism (see, e.g., |17j for a proof). 

Let us fix some arbitrary q (e.g., see the inset). The first observation is 
that the "corner" triangle, namely the triangle touching the origin does 
not flip orientation under x q - This can be seen by using the fact that \"f~\0 < j. This actually implies that 
all triangles with an edge on the segments [0, £] and [0, rf\ (colored red in inset) are not flipped. 
We are left with the triangles of the type colored green and gray in the inset. Since, as we noted before, 
scaling does not change the orientation we can consider w.l.o.g. triangles of the form A(l, 1 + 77, 1 + £) £ T q . 
Since every triangle in T q is similar to T, we can further assume that the corner at 1 corresponds to in T. 
Since 6 < 60.4° < n/2 one must have Re(£)Re(?7) > 0. The triangles for which Re(£) < and Re(^) < 
(corresponds to the gray triangles) are not flipped since h(l) = 1, h(l + £) = (1 + £) 7 , |(1 + £) 7 | < 1, and 
arg((l + 7 ) < f since 6\j\ < f and arg(l + £) < 9. Similarly, arg((l + 77) 7 ) > -§ and therefore the 
triangle A(l, h(l + rj), h(l + £)) has the same orientation as A(l, 1 + r), 1 + £). 

We are left with the triangles for which Re(£) > and Re(r/) > (the green triangles in the inset). We will 
later use the fact that since all the triangles are isoceles a consequence from the cosine law is that, for all 
"green" triangles, 

(A ' 1} K 1 - /onl m > 

y/2(l + cos(6>)) 

and similar bound holds for \rj\. 

A sufficient condition that will prevent triangle A(l, 1 + 77, 1 + £) from flipping is 

Re(l + £) 7 - 1 > , Re(l + t?) 7 - 1 > 0. 
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We will show the first inequality, and the second is proved similarly. In the following we use the binomial 
expansion 2J. Denote £ = re lip , then 

R*(l+0 7 -l = 7Rce+ 7( VW + 7(7 ~o, (7 ~ 2) IK 3 1 



7r cos(iy5) 



3! 

lh^l r 2 cos(2 ^) + 7(7-l)(7-2) r3 cQs{3ip) 



7 — 1 

jr ( coa(ip) H — rcos(2<^) + 



3! 

(7~l)(7-2) 
3! 



r 2 cos (3^) 



(A.2) 



where dj = ^ 7 ^.3? 1 f° r J — 1, 2, and do = 1. Let k G N be such that fc < 7 < fc + 1 . Since 

0<7 — fc<lwe have 



\a k 



( 7 -l)(7-2)-...-( 7 -fc) 



2-3-...-(fc + l) 
Using induction one proves that for j = k + 1, k + 2, ... 

1 

1 " " 3+1 



< 



k + 1 



We can now bound (A.2 1 as follows 



Re(l + £) 7 - 1 = jr aj cos ((j + l)(p)r j + o 3 - cos ((j + l)ip) r j 



> 



> 



7 



(cosM- ~^rA 
\ 3 =k+i J ^ J 



> 7r ( cos(y) (— log (1 — r) — r) 



> ir I cos(0) + ^1^1) + 1 



(A.3) 

where in the second row we used the fact that for 1 < j < k, j + 1 < [7], and therefore cos((j + 1)^) > 
cos(|y|6>) > as 6Tt1 < ^/2. Plugging in ( |A.3| we get 

log(l-[2(l + cos(#))r 1/2 ) ^ 



Re(l + £) 7 -l > 7 r ^cos(<9) 
= irT{9), 



[2(l + cos(0))] 



-1/2 



Investigating the function T(8) (see inset for its graph, in red) 
reveals that for 9 < 60.4° T(0) > and therefore Re(l + 7 - 
1 > 0. Similarly, Re(l + if) 1 — 1 > and therefore the triangle 
A(l, h(\ + r]),h(l + £)) is not flipped. □ 
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